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Abstract 



Predictions of two channels in the three-body decays of the charmed mesons are made 
within the heavy hadron chiral perturbation theory. There still exists the problem that 
the theoretical expectation is too small compared to the experimental data. 
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Nonleptonic weak decays of the charmed mesons have been studied extentively 
in the past two decades. Previous theoretical studies are focused on the two-body 
cases [0]. The experimental measurements have also been achieved in some three-body 
channels ||. The experimental results, however, are not well understood quantatively 
due to two reasons. On the theoretical side, there exists no method in the literature 
which allows one to perform the calculation reliably. On the experimental side, there 
are so many open resonant channels which contribute to the final three-body states 
that the data are difficult to be analysed. 

In the present work, we use the heavy hadron chiral perturbation theory|| |], 
Q(HCHPT) to study the non-resonant three-kaon decays of the charmed mesons. In 
fact, application of chiral perturbation theory in this kind of study is not a new idea. 
In the past, the U(A)i ® U(4)r chiral symmetry has been used in |], [jj. Because this 
symmetry is badly broken, these predictions are totally not under control. 

HCHPT introduced in || can be described in the following. The QCD lagrangian 
for the light quark (u, d, s) sector possesses the SU(3)l <8> SU(3)r chiral symmetry. 
The heavy quark (b or c), which transforms as singlet under the chiral symmtery, has 
the spin-flavor symmetry in the limit that its mass is taken to be infinity [§]. As the 
consequence, the two lowest lying mesons containing one heavy quark are degenerate 
in the heavy quark limit, and can be expressed by the superfield 

H a = ^{l+/b){D^f-D alB ), (1) 

where use of the charmed mesons, cq a with a = 1,2,3 corresponding to D 0( -*\ D +( -*\ 
D+W, has been made as the example. In ([[]) we have suppressed the explicit depen- 
dence of H a on its velocity v. The strong interactions of the heavy mesons with the 
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pseudo-Goldstone bosons tt, K and 77 at low energy can be constructed by taking the 
derivative expansions on the pseudo-Goldstone field S = £ 2 = exp(2iM/ f), where 
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and / is the decay constant of pseudo-Goldstone bosons. In the derivative expansions, 
higher order terms are suppressed by powers of 1/Acsb with the chiral symmetry break- 
ing scale Ac sb ~ 1.2GeV from the naive dimensional analysis 0. As the superfield 
(H) is used, the requirements of the heavy quark symmetry is satisfied automatically. 
Higher order terms which violate the heavy quark symmetry are suppressed by powers 
of 1/Mq and can be incorporated into HCHPT [[lO| . 

To the leading order in both the derivative and the 1/Mq expansions, the effective 
lagrangian in HCHPT is 



C = -iTrH a v^H a + \irrH a H h v^ + d^ + t,d^ 
+\igTrH a H bltll s(t+d^ - $^+) h + • • • , 



(3) 



where the trace is taken over the spinor space. The coupling g in (|]) is estimated to be 
of order one and can be extracted from the partial width of the strong decays D* — ► Dtt. 



Up to now it has only an upper bound \g\ < 0.63|11[. We will use g = 0.5 ± 0.1 in the 
numerical estimations. 



In HCHPT, semiloptonic decays of heavy-to-light transitions are descibed by the 
effective weak current 



(4) 



where both sides transform as (3^, 1r) under the SU(3)l <8> SU(3)r chiral symmetry, 
and 

a = fD\frriD- (5) 
The light quark current is described in the same way as in the usual chiral perturbation 
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ft7 /i (l-7te)(& = E ! f^ S * Et W ( 6 ) 



The effective hamiltonian for Cabibbo-allowed three- K decays of the charmed mesons 



is: 



G 

n eff = -j=V cs V: d [ ai u^(l - 75)^(1 - 75 )c + a 2 s 7 ^(l - 75)^(1 - 75 )c], (7) 

where a\ = 1.2 and a 2 = —0.5, which are the most favored values in the phe- 
nomenological analyses in the two-body decays[|l3|, will be used in numerically es- 
timations. In dealing with the nonleptonic decay ampitudes we use the factorization 



ansatz[]13] under which the amplitudes for the three-body decays depicted in Figure 1 
for D° -> K+K-R and D+ -> K + K°K° are: 

<K+( q+ )K-(q^K% qo )\H eff \D > = ^V cs V: d 

[ax < K+{q + )K Q {q Q )\u^{l- lb )d\Q >< K~ {q^s^l - 75 )c|£>° > 
+a 2 < K°(q )\sr(l - l5 )d\0 >< K+(q + )K-(q_)\u lfl (l - <y B )c\D° >] 



<K+{ q+ )K\ qi )K%q 2 )\H eff \D+> = ~^V cs V: d] 



G_f 

^2' cs ' ud y 2 

[ai < K + ( q+ )K°( qi )\ur(l-ls)d\0 >< K°(q 2 )\s^(l-j 5 )c\D+ > 
+a 2 < K°( qi )\s^(l - j 5 )d\0 >< K + ( q+ )K ( q2 )\uY(l - l 5 )c\D+ > 
+ qi Q2 ]• 

In Figure 1 we have discarded the W-exchange and the W- annihilation diagrams which 
are expected to be highly suppressed. 



The applicability of the chiral perturbation theory in D — > 3K lies in the following 
reasons. In the final states, the maximum energy of each of the .fT-meson in the rest- 
frame of the D meson is 

E max ~ 0.73GeV, (9) 
while the maximum value of the invariant mass of any two i^-mesons is 

{^Mf^max = m D -m K ~ 1.3GeV, (10) 

which is a little larger than the estimation of Acsb ~ 1.2GeV from the naive dimen- 
sional analysis0. However, Acsb can be also taken as 1.5GeV, as has been analysed in 
the literature [Q. In this case, the whole phase space of these decays are in the region 
where HCHPT is applicable. On the other hand, even if Acsb ~ 1.2GeV is taken, the 
phase space where HCHPT can be applied is still dominant, because it corresponds to 
the much large area in the Dalitz plot. Note that discarding of the annihilation dia- 
gram is important to avoid the terms proportional to the invariant mass of the three 
final particles. 

Note that the two channels depicted in Figure 1 are the only three-body ones which 
can be analysed in HCHPT. The corresponding hadronic matrix elemnets in (§) are 
estimated by calculating the Feynman diagrams in HCHPT which are depicted in 
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Figure 2. The results are 

<K+(q + )K\q )\uY(l-l5)d\®> = (?+ ~ *>)", 

<A-°(g )|^(l-75)d|0> = i<tf/, 

<K-(g_)|s 7M (l-75)c| J D°> = Y 1 (q_), + Y 2 (q_)^, 

<K + (g + )^-(g_)|« 7 M(l- 75)c|£>° > = X 1 (q+, + X 2 (g+, 

+ ^3(g+,g-) M + x 4(g+>^)M' 

<ir + (g + )K°( gi )|M7 M (l-75)rf|0> = (<z+-<zi)", 

<K ( gi )|s7 M (l -75^0 > = </, 

<^°(g 2 )|s7 M (l-7 5 )c| J D + > = Y 1 (q 2 ), + Y 2 (q 2 )„ 

<K + (q + )K°(q 2 )\u llI (l- l5 )c\D+> = X 1 (q + , g 2 ) M + X 2 (q + , g 2 ) M 

+ ^s(g+,g2) M + X 4 (g + ,g 2 ) M , 



where 



v / \ JdPd,* v- (q-q+) 

Xi{q+,q)t, = 1- 



P 2v(q + q + Y 



~y / \ ■ 2 /o P -Dm q-q+-v-qv-q+ 

X 2 (q+,q)^ = -%g A 



X^iq+.q)^ = g- 



P v ■ (q + q+)(v ■ q + A)' 
ifD s -rn D q^ + v ■ qP D ^ 



P v ■ q + A 



X4{q+,q)^ = 2 j 2 
coming from Figure 2(a)-(d), respectively, and 

v(n \ _ Jds -mog^ + v ■ qP Dli 

Y 2 (q), = ~Y Pd » 
from Figure 2(e)-(f). We have denoted 

A = m D * — mr>- 
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In the numerically evaluations, we take fu = fu s = 0.2GeV and f = fx — 
0.161GeV. The effective coupling g is taken to be 0.4, 0.5 or 0.6(the corresponding 
formfactor at the maximum momentum transfer in the semileptonic decays of D — ► K 
is \ f+(q 2 m )\ = 1.08, 1.20 or 1.38, while the experimental value is 1.30±0.5 if a monopole 
behavior of the g 2 -dependence is used|L5|]). We give our results in Table 1, together 



with the comparisions with both the estimations from U(4)l <8> i7(4)#J7| and the ex- 
perimental data||15|| . Note that no numerical prediction has been made in for the 



two channels we have studied. As has been found in the U(A)i ® U(4) R studies]?]], 
there are some three-body channels whose measured branching ratios are larger than 
the theoretical expectations by more than one order. In the two channels studied in 
the present work, we are still suffered from the same problem even if our calculations 
are based on more reliable foundation. This problem cannot be solved by going to the 
higher order expansions in HCHPT, otherwise the expansions will not converge. It is 
also impossible to attribute this problem to the omissions of the W-annihilation and 
the W-exchange diagrams because they are suppressed compared to those in Figure 1. 

To bridge the gaps between the theoretical estimations and the experimental mea- 
surements, further studies at the future r — Charm factory are essential, where strict 
subtractions off the contributions from many resonant channels should be carried out. 
In the meantime, the interference effects between resonant and non-resonant channels 
are also needed to be studied by both the theoriests and the experimentists. 

The author would like to acknowledge G. Eilam for suggestion of the present work 
and helpful discussions. This research is supported in part by Grant 5421-3-96 from 
the Ministry of Science and the Arts of Israel. 
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Table 



Table 1 Comparisions of numerical results. Our results using different values of g are 
given in the second column. 



process 


g =0.4 g=0.5 g=0.6 


M 


Exper. 15 




Br(D° -> K + K~K°) 


1.3 x 10~ 4 1.8 x 10~ 4 2.3 x 10~ 4 


6.3 x icr 5 


(4.9 ±0.9) x 10^ 3 
(non — 0) 


Bi(D + -> K+K°K°) 


6.4 x lCT 4 8.3 x io- 4 1.1 x 1CT 3 




(3.1 ±0.7)% 



Figures 



Figure 1 The Feynman diagrams for the D° — > K° and D + — > K + K°K°. 



Figure 2 The Feynman diagrams used in HCHPT to calculate the hadronic matrix 
elments between the heavy and the light mesons. 
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